Mott transitions in the Hubbard model with spatially-modulated interactions 



Akihisa Koga and Takamitsu Saitou 
Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 

Atsushi Yamamoto 
RIKEN, Advanced Institute for Computational Science, 7-1-26, 
Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan 
(Dated: July 11, 2012) 

We study two-component fermions in optical lattices with spatially alternating on-site interactions 
by means of dynamical mean-field theory. Calculating a quasiparticle weight, double occupancy, and 
order parameters for each sublattice, we discuss low-temperature properties in the system. When 
the repulsive interactions alternate, the magnetically ordered state is realized at half filling, while 
the superfluid state is realized in the attractive case. Furthermore, we find that the metal-insulator 
transition occurs against magnetic and superfluid fluctuations when the repulsive and attractive 
interactions alternate. In the doped system, the commensurability emerges due to the repulsive 
interactions and the insulating state is realized down to the quarter-filling. 



I. INTRODUCTION 

Recent extensive experimental and theoretical inves- 
tigations on ultracold gases have providing a variety of 
interesting topics. Typical examples are ultracold atoms 
loaded on optical lattices with tunable and controllable 
parameters^ where many remarkable phenomena have 
been observed such as the Mott transitions^— and BCS- 
BEC crossover £~— Recently, the spatial modulation of the 
atom-atom contact interaction has been realized in 174 Yb 
gas system^ which stimulates further theoretical inves- 
tigations on particle correlations in the ultracold atomic 
systems. 

One of the interesting questions is how the spatial 
modulation in the interactions affects low-temperature 
properties in fermionic optical lattice systems. Here, we 
consider the optical lattice with alternating on-site in- 
teractions as a simplest model. This should capture the 
essence of the systems with the modulated interactions 
with a finite period. In our previous study,— we have 
clarified that in the half-filled system, the Mott transi- 
tions occur when the magnitudes of two interactions in- 
crease. However, it is naively expected that magnetically 
ordered state or superfluid state is stabilized in the sys- 
tem only with repulsive or attractive interactions. There- 
fore, it is necessary to discuss whether or not the Mott 
transition is indeed realized against magnetic and super- 
fluid fluctuations. It is also interesting to clarify how the 
alternating interactions affects low-temperature proper- 
ties in the doped system. This should be important in 
discussing low-temperature properties in realistic opti- 
cal lattice systems where the local particle density varies 
due to the existence of the harmonic potential)^ ' 10 ' 11 
Moreover, the system with alternating interactions can 
be regarded as the natural extension of the two-band 
modclfi^— which has been studied in condensed mat- 
ter physics. Therefore, it is instructive to systematically 
study how low-temperature properties are affected by the 
spatial modulation of the on-site interactions. 

In this paper, we consider the Hubbard model with al- 



ternating on-site interactions to discuss how particle cor- 
relations affect low-temperature properties. To this end, 
we make use dynamical mean-field theory (DMFT}i£~— 
with a continuous-time quantum Monte Carlo (CTQMC) 
method^ This combined method has an advantage in 
dealing with magnetically ordered and superfluid states 
on an equal footing in the strong coupling region. Thus, 
we will show how such states are affected by alternating 
interactions. Furthermore, we discuss the stability of the 
metal-insulator transition in the system. 

Our paper is organized as follows. In Sec. [Hi we intro- 
duce the model Hamiltonian and briefly summarize our 
theoretical approach. In Sec. IIII1 we demonstrate how 
the magnetically ordered and superfluid states are real- 
ized at half filling. It is clarified that the metal-insulator 
transition occurs in the system with both repulsive and 
attractive interactions. We also discuss the effect of the 
hole doping to clarify the role of the repulsive interac- 
tions in Sec. IIVI A brief summary is given in the last 
section. 



II. MODEL AND METHOD 

We consider two-component fermions in an optical lat- 
tice with alternating on-site interactions, which should 
be described by the following Hubbard Hamiltonian on a 
bipartite lattice, 

H = -t c L c j<y ~ + ha) n i° 



2j' 



(1) 



where cj^Cio-) creates (annihilates) a fermion at the «th 
site with spin a and = c\ a Ci a . t is the hopping in- 
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tegral and U a is the on-site interaction at the sublat- 
tice a(= A, B). \i is the chemical potential and h is 
the magnetic field. We note that in our model Hamil- 
tonian cq. ([1}, (n) = 1 (half filling) is realized when 
/i = h = 0, where n = J2ia n i<?/N an d N is the to- 
tal number of sites. In the case, there exists the sym- 
metry, in addition to the trivial symmetry axis (Ua = 
Ub)- By applying the particle-hole transformation^ 
Cj-f — > Cj-j- and — > (— l) 1 ^ to the model Hamilto- 
nian with (Ua,Ub, H,h), we obtain the model Hamil- 
tonian with (—UA,—UB,h,fJ,). Therefore, in the case 
with fj, = h = 0, the Hubbard models with (Ua,Ub), 
(Ub,Ua), {—Ua,~Ub) and (—Ub,—Ua) are identical. 
In the repulsive model (Ua > and Ub > 0), a mag- 
netically ordered state expected naively is equivalent to 
the superfluid state in the attractive model (Ua < and 
Ub < 0). We also note that the particle number imbal- 
ance (ua ^ tib); which may be regarded as the density 
wave state, generally appears since the system has a two- 
sublattice structure, where n a = 2 J2iea a n ia/N. 

The Hubbard model with alternating interactions has 
been studied in one-dimensional systems, where the pos- 
sibility of anomalous metallic states is discussed^ In 
the paper, we deal with the infinite-dimensional Hubbard 
model to discuss how alternating interactions affect low- 
temperature properties. To clarify this issue, we employ 
the DMFT method^i 9 . In DMFT, the lattice model is 
mapped to an effective impurity model, where local parti- 
cle correlations are taken into account precisely. The lat- 
tice Green's function is then obtained via self-consistent 
conditions imposed on the impurity problem. The treat- 
ment is formally exact in infinite dimensions, and even in 
three dimensions, DMFT has successfully explained in- 
teresting physics such as the Mott metal-insulator tran- 
sition. 

When DMFT is applied to the system with a sublatticc 
structure, the Green's function is given as^ 

G(k,z)- 1 = G (k,z)- 1 -t(z), (2) 

with 

G (k, z) = ( {Z + Q* + » A *> , ^ ~ £k& l ) 1 , 

V ; V -£kCT 2 (z + h)<T + HBVz J 

(3) 

and 

fi W-( fiA o W fi. (,))- « 4 » 

where &q is the identity matrix, a z is the z-component 
of the Pauli matrix, [i a = fx — U a /2 and £k is the dis- 
persion relation for the bare band. S Q (z) is the self- 
energy for the ath sublatticc in the Nambu formal- 
ism. The local lattice Green's function is obtained as, 
G a (z) = J dkG aa (k, z). Here, we use a semicircular den- 
sity of states, p(x) = -£p\J 1 — {%) 2 where D is the half- 
bandwidth, which corresponds to an infinite coordination 



Bethe lattice. The self-consistency equation is then given 
by 

Q a (z) = za + fi a a z - f^j a z Ga(z)a z , (5) 

where Q a (z) is the non-interacting Green's function of the 
effective Anderson impurity model for the ath sublatticc. 

There are various numerical methods to solve the ef- 
fective impurity problem such as the iterative perturba- 
tion theor y 24 ' 25 and the exact diagonalizationJ^ Here, 
to discuss low-temperature properties quantitatively, we 
make use of the CTQMC method? 2 ^, This technique 
has recently been developed^ and has successfully been 
applied to general classes of models such as Hubbard 
modelr^— periodic Anderson model) 29 ' 30 Kondo lat- 
tice model 3 -^ and Holstein-Hubbard model^ Here, we 
use the hybridization-expansion version of the CTQMC 
method^ extended to the Nambu formalism^ 3 - This al- 
lows us to directly access the superfluid state at low tem- 
peratures. In our CTQMC simulations, we measure nor- 
mal and anomalous Green's functions on a grid of one 
thousand points. Furthermore, we also use the numeri- 
cal renormalization group (NRG) method 3 ^ to discuss the 
stability of the normal metallic state complementary. In 
the NRG, one discretizes the effective bath on a logarith- 
mic mesh by introducing a discretization parameter A. 
The resulting discrete system can be mapped to a semi- 
infinite chain with exponentially decreasing couplings, 
which allows us to access and discuss properties involv- 
ing exponentially small energy scales^ To ensure that 
sum rules for dynamical quantities are fulfilled, we use 
the complete-basis-set algorithm^ - — We observe that 
obeying the sum rules is mandatory to properly describe 
low-energy properties of the system away from half fill- 
ing. In the NRG calculations, we use a discretization 
parameter A = 2 and keep four thousands states at each 
step. 

In this paper, we use the half-bandwidth D as the unit 
of energy. To discuss low-energy properties of the sys- 
tem with alternating interactions systematically, we cal- 
culate quasiparticle weight z a = (1 — dRcT, a (uj) / duj)^ 1 , 
double occupancy D a = (n a ^n a ^), magnetization m a = 
(n-at — n a i)/2, and pair potential A a = (c a ^c a {) for 
each sublatticc. When the CTQMC method is used 
as an impurity solver, we calculate the quantity z a = 
(1 — ImEo,(iwo)/wo) _1 as a quasiparticle weight at finite 
temperatures, where ujq = ir/f3. Furthermore, by ap- 
plying the maximum entropy method 3 ^— to the Green's 
functions, wc deduce the spectral functions, which allows 
us to discuss static and dynamical properties of the sys- 
tems. 



III. RESULTS AT HALF FILLING 

We first discuss low-tcmperaturc properties in the 
Hubbard model at half filling (n = When the signs 
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of two interactions arc the same (UaUb > 0), the or- 
dered ground-state is naively expected at zero tempera- 
ture. Namely, the magnetically ordered state is stabilized 
in the repulsive case, while the superfluid state is in the 
attractive case. Here, we consider the Hubbard model 
with repulsive interactions since the attractive model at 
half-filling is equivalent under the particle-hole transfor- 
mation, as discussed above. To clarify how the magneti- 
cally ordered state is affected by alternating on-site inter- 
actions at low temperatures, we calculate the staggered 
magnetization ttlaf = ^T m i/^- The obtained re- 

sults at T/D = 0.05 are shown in Fig. [TJ In the weak 
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FIG. 1: (Color online) Density plot of the spontaneous mag- 
netization in the system with T/D = 0.05. 



coupling region, the paramagnetic metallic state is real- 
ized with ttiaf = 0. Increasing the interactions at both 
sublattices, the staggered magnetization is induced and 
the phase transition occurs to the magnetically ordered 
state. In the state, the magnitude of magnetization is 
almost uniform in the system even with alternating in- 
teractions (Ua 7^ Ub) although the Green's functions 
strongly depend on the sublattices. Further increase in 
the interactions leads to another phase transition to the 
paramagnetic state since thermal fluctuations are rela- 
tively enhanced in comparison with magnetic fluctua- 
tions. Since the effective intersite coupling should be 
scaled as t 2 / (Ua + Ub), the phase boundary is almost lin- 
ear (Ua + Ub = const.) at finite temperatures, as shown 
in Fig. Q] It is expected that the magnetically ordered 
region is extended at lower temperatures. 

On the other hand, when the signs of interactions are 
different from each other (UaUb < 0), ground state 
properties should be nontrivial. In this case, a non- 
magnetic ground state is expected since the cooperative 
phenomenon is hard to be realized in the framework of 
DMFT. In fact, we could not find any signature of the 
ordered state down to T/D = 0.02. To discuss the stabil- 
ity of the normal metallic state, we calculate the quasi- 
particlc weights for both sublattices with a fixed ratio 



Ub/Ua = -2 at T/D = 0.02, as shown in Fig. [2] When 




FIG. 2: (Color online) Quasiparticle weight for each sublat- 
tice as a function of Ua/D with a fixed ratio Ub/Ua = —2. 
Solid lines with circles are obtained from DMFT with the 
CTQMC method at T/D = 0.02. Dashed lines with squares 
are obtained from DMFT with the NRG method. Inset shows 
the double occupancy in the system. 

the interactions are turned on, the quasiparticle weights 
za and zb decrease from unity in slightly different way 
reflecting the difference of on-site interactions. Namely, a 
strong renormalization appears in the sublatticc B. Fur- 
ther increase in the interactions induces the jump sin- 
gularity in both curves around Ua/D ~ 1.7. Then the 
double occupancy approaches zero (half) at the sublat- 
tice A(B), as shown in the inset of Fig. [5] These imply 
that the first-order phase transition occurs to the insulat- 
ing state, where the singly occupied states with spin are 
realized at the sublattice A, while the empty or doubly 
occupied states are equally realized at the other. There- 
fore, we conclude that the Mott and pairing transitions 
simultaneously occur in the sublattices A and B, respec- 
tively. 

The density of states at each sublattice at T/D = 0.02 
is shown in Fig. [3] It is found that two kinds of broad 
peaks appear in the high energy region. One of them 
is located around w ~ ±£/ Q /2, which forms the Hub- 
bard gap due to the local Coulomb interactions. The 
other broad peak around lj ~ ±C/ a /2 is induced by the 
on-site interaction at nearest neighbor sites although it 
may be invisible on this scale for the density of states 
Pa- When Ua/D = 1.75, the sharp quasiparticle peaks 
for both sublattices appear in the vicinity of the Fermi 
level, which implies that the system is in the metallic 
state close to the Mott transition point. On the other 
hand, when Ua/D = 2.0, the gap structure appears in 
both sublattices, which means that the insulating state is 
induced by the Mott and pairing transitions. Therefore, 
we conclude that the single metal-insulator transition in- 
deed occurs in the system with alternating interactions 
(UaU b < 0). 

To examine the nature of the phase transition at zero 
temperature, we also use the NRG method as an impu- 
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contrast, when the signs of the interactions are different 
(V 'aUb < 0), any ordered states should not be stabilized 
and the metal-insulator transition occurs. We find that 
the obtained critical points are in good agreement with 
the phase boundary U aUb = — (3-D) 2 obtained from the 
linearized DMFT method.— The Mott insulating state is 
stabilized in the strong coupling region, where the singly 
occupied states are realized in the sublattice with the re- 
pulsive interaction, and the empty or doubly occupied 
states in the other. In the following, we discuss how al- 
ternating interactions affects low-temperature properties 
in the doped system. 

IV. HOLE DOPING EFFECTS 



FIG. 3: Solid (dashed) lines represent the density of states for 
the sublattice A(B). Upper (lower) panel shows the results for 
the metallic state with Ua/D — 1.75 (the insulating state with 
Ua/D = 2.0) in the system with U B /Ua = -2 at T/D = 0.02. 



rity solver. In the weak coupling region, the NRG re- 
sults are in good agreement with the CTQMC results, 
as shown in Fig. [2] However, increasing the interac- 
tions, the quasiparticle weights do not have the jump 
singularity, and smoothly reach zero at the same critical 
point (Ua/D <~ 2.1). This suggests the existence of the 
second-order quantum phase transition in the half-filled 
Hubbard model with alternating interactions, which is 
similar to the nature of the conventional Mott transition 
in the paramagnetic state. 

By performing similar calculations, we obtain the zero 
temperature phase diagram, as shown in Fig. |4j When 
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FIG. 4: Ground-state phase diagram of the half-filled Hub- 
bard model with alternating interactions. Solid lines repre- 
sent the phase boundaries obtained from the linearized DMFT 
method^ and circles are obtained from DMFT with the NRG 
method. 



Ua and Ub are positive (negative), the magnetically or- 
dered (superfluid) ground state should be realized. In 
this case, the Mott (pairing) transition never occurs. By 



In the section, we discuss a hole doping effect on the 
system with alternating interactions. At first, we con- 
sider the system only with repulsive interactions. It is 
known that the phase separated magnetically ordered 
state is realized in the conventional Hubbard model close 
to half filling at zero temperature^ To study how the 
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FIG. 5: (Color online) Density plot of the spontaneous mag- 
netization (a) and the particle density in the sublattice A (b) 
for the doped system (n = 0.95) at T/D = 0.05. 



spatial modulation in the interactions affects the mag- 
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netic structure at low temperatures, we calculate the 
magnetization tuaf m the system with n = 0.95, as 
shown in Fig. [5] (a). The magnetization is monotoni- 
cally decreased by the hole doping, in comparison with 
the half-filled case (see Fig. [TJ. It is also found that the 
phase boundary in the strong coupling region somewhat 
curves. This behavior should be related to the particle 
occupation in each sublattice. Figure [5] (b) shows the 
particle number for the sublattice A. It is found that the 
commensurability emerges in the sublattice A (tia ~ 1) 
when Ua ^ Ub- I n the region, magnetic correlations 
are enhanced easily, which results in the curved phase 
boundary. Further hole doping suppresses magnetic cor- 
relations, which leads to the phase transition to the para- 
magnetic stated 

By contrast, it is known that the supcrfluid state is 
stable against the hole doping. To clarify how the al- 
ternating interactions affect the superfluid state, we deal 
with the attractive Hubbard model with n = 0.75 and 
Ub/Ua = 5. Calculating the particle density and pair 



U B /U A =5 




FIG. 6: (Color online) Particle density for each sublattice (a) 
and pair potential (b) as a function of Ua/D with a fixed 
ratio Ub/Ua = 5. Solid (open) symbols are the results for 
sublattice A(B) obtained from DMFT with CTQMC method. 



potential, we obtain the results at fixed temperatures 
T/D = 0.02 and 0.05, as shown in Fig. [BJ The para- 
magnetic metallic state is realized below a certain crit- 
ical interaction (Ua) c , and the particle number imbal- 
ance appears with ua ^ rig, as shown in Fig. [S] (a). 
In this case, the Hubbard gap (<~ \Ub\) appears in the 
density of states for the sublattice B, while a low-energy 



peak clearly appears in the other, as shown in Fig. [7] 
Beyond the critical interaction (Ua)c, the pair potential 
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FIG. 7: (Color online) Density of states when Ua/D = 
-0.4, -0.5 and -0.75 at T/D = 0.05. 

A is induced and the superfluid state is realized. Then 
the superfluid gap appears in the vicinity of the Fermi 
level, as shown in Fig. [7] By examining critical behav- 
ior A ~ \Ua — (Ua) c \ with the exponent j3 = 1/2, we 
determine the critical interactions (Ua/D) c = 0.47 at 
T/D = 0.05 and (U A /D) C = 0.27 at T/D = 0.02, as 
shown in the inset of Fig. [5] (b). It is also found that the 
particle number imbalance is suppressed once the system 
enters the superfluid state, as shown in Fig. [HI (a). Since 
this tendency is clearly found at lower tempertures, we 
can say that the realization of the superfluid state makes 
the system uniform. Similar behavior can be found in 
the half-filled Hubbard model with a uniform attractive 
interaction where the hole doping destabilizes the density 
wave state and the genuine superfluid state is stabilized. 

When UaUb < 0, interesting behavior appears in the 
doped system (n = 0.75). The results for the system with 
a fixed ratio Ub/Ua = —2 at T/D = 0.02 are shown in 
Fig. [S] The increase in the interactions monotonically 
decreases the quasiparticle weight for each sublattice and 
a singularity appears at Ua/D ~ 1.7, which is similar to 
the half-filled case. However, a further increase in the 
interaction leads to different behavior. In the sublattice 
A, the quasiparticle weight becomes almost zero and the 
local particle density is close to the half filling, as shown 
in the inset of Fig. [51 Our results suggest the existence of 
the Mott transition in the sublattice A. In fact, the Mott 
gap clearly appears in the density of states, as shown 
in Fig. [9l By contrast, in the sublattice B, the local 
particle density becomes away from the half filling. Then 
the quasiparticle weight increases beyond the transition 
point and the gap structure appears around the Fermi 
level. This behavior is characteristic of the spin-gap state 
induced by the pairing transition in the doped system.— 
Therefore, we can say that even in the doped system, the 
Mott and pairing transitions occur at finite temperatures. 
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FIG. 8: (Color online) Quasiparticle weight for each sublat- 
tice as a function of Ua/D with a fixed ratio Ub/Ua — —2. 
Solid lines with circles are obtained from DMFT with the 
CTQMC method at T/D = 0.02. Dashed lines with squares 
are obtained from DMFT with the NRG method. Inset shows 
the particle density for each sublattice. 
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FIG. 9: (Color online) Density of states for the system with 
U B IV a = -2 at T/D = 0.02 when U A /D = 1.6 (upper panel) 
and 1.87? (lower panel). 



However, the nature of the phase transition at zero 
temperature may not be trivial since in the system with 
the uniform interactions, the Mott transition is of second 
order at half filling and the pairing transition is of first 
order away from the half filling. To clarify the stability 
of the normal metallic state, we calculate the renormal- 
ization factors at zero temperature by means of the NRG 
method. We find that the quasiparticle weights smoothly 
decrease up to Ua/D ~ 2.0. In the case, a large difference 
appears between the NRG and CTQMC results, contrast 
to the half-filled case. This originates from the fact that 
the particle number at each sublattice is sensitive to the 
interaction strength and thermal fluctuations. When the 
interaction is large enough, oscillation behavior appears 
in the DMFT iterations and wc could not find the stable 



solution. This may imply that the first-order phase tran- 
sition occurs and the phase separation is realized at zero 
temperature, which has been discussed in the attractive 
Hubbard model4^ This point thus remains an open issue 
for future studies. 

To clarify the role of the attractive and repulsive inter- 
actions in the doped system, we also consider the system 




FIG. 10: (Color online) Quasiparticle weight for each sublat- 
tice as a function of Ua/D with a fixed ratio Ub/Ua = —1/2. 
Solid lines with circles are obtained from DMFT with the 
CTQMC method at T/D = 0.02. Dashed lines with squares 
are obtained from DMFT with the NRG method. Inset shows 
the particle density for each sublattice. 

with Ub/Ua = —1/2, where the magnitude of the at- 
tractive interactions is smaller than the other. Fig. [TU1 
shows that when Ua/D ~ 1, the renormalization fac- 
tor in the sublattice A is rapidly decreased, while the 
other is little decreased. This behavior is contrast to 
the case Ub/Ua = —2 (see Fig. [5]), which means the 
absence of the symmetry axis (Ua = —Ub)- On the 
other hand, when the system approaches the phase tran- 
sition point, the particle number for the sublattice A with 
repulsive interactions approaches half and both quasi- 
particle weights decrease. At last, the jump singular- 
ity appears in the curves of the quasiparticle weights at 
Ua/D ~ 3.3, where the Mott and pairing transitions oc- 
cur simultaneously. Therefore, the nature of the phase 
transition in the doped system is essentially the same as 
the half-filled case. We wish to note that this is contrast 
to low-temperature behavior in the doped system with 
both repulsive interactions. In the system, the stronger 
repulsive interaction tends to recover the commensurabil- 
ity in the corresponding sublattice, but the weaker one 
never induces the phase transition in the intermediate 
band filling. Then the heavy metallic state is realized in 
the strong coupling regime. Therefore, we can say that a 
single phase transition discussed here is stabilized by two 
effects: the Mott transition together with the commensu- 
rability due to the repulsive interaction and the pairing 
transition induced by the attractive interaction. 

When the particle number is smaller than the quarter 
filling (0 < n < 0.5), it is expected that the common- 



7 



surability never emerges and no phase transition occurs. 
To confirm this, we show the renormalization factors in 
the systems with n = 0.4 and 0.6 at T/D = 0.02, as 
shown in Fig. QTJ In the case n = 0.6, the singularity 




FIG. 11: (Color online) Quasiparticle weight for each sublat- 
tice as a function of Ua/D with a fixed ratio Ub/Ua = —2 
at T/D = 0.02. Circles (triangles) are the results for n — 
0.6(0.4), obtained from DMFT with the CTQMC method. 



appears around Ua/D ~ 1.6, which suggests the exis- 
tence of the single phase transition. On the other hand, 
when n = 0.4, no singularity appears in the curves of the 
renormalization factors. Therefore, in the doped system 
with n < 0.5, the commensurability never emerges and 
no Mott transition occurs. 

In this paper, we have discussed how the spatially- 
modulated interactions affect low-temperature proper- 
ties in the optical lattice system. We have determined 
the rich phase diagram of the system. Although this 
model may be one of the realizable toy models, we have 
clarified the role of the repulsive and attractive inter- 
actions. In particular, the repulsive interaction plays a 
crucial role in stabilizing the metal-insulator transition 



in the doped system. This is similar to the essence of the 
orbital-selective Mott transition in the doped multiband 
system , 45 i 46 where the commensurability in the narrower 
band emerges due to the repulsive interactions. It is an 
interesting problem whether or not such phase transition 
indeed occurs in the finite dimensions, which is now un- 
der consideration. 



V. SUMMARY 

We have investigated the fermionic optical lattice sys- 
tem, which is described by the Hubbard model with al- 
ternating on-site interactions. By combining dynamical 
mean-field theory with the continuous-time Monte Carlo 
method, we have studied low-temperature properties sys- 
tematically, calculating the quasiparticle weight, double 
occupancy, and order parameters for each sublattice. It 
has been clarified that the magnetically ordered state is 
realized in the system with only repulsive interactions, 
while the superfluid state is realized in the attractive 
case. We have also found that when the repulsive and at- 
tractive interactions alternate, the metal-insulator transi- 
tions occur against magnetic and superfluid fluctuations. 
In the doped system, we have found that the commen- 
surability emerges due to the repulsive interactions and 
the insulating state remains with n > 1/2. 
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